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The Fermi Bubbles have been imaged in sub-TeV gamma rays at Fermi-LAT, and, if their origin 
is hadronic, they might have been seen with low statistics in ~ 0.1 — 1 PeV neutrinos at IceCube. 

We discuss the detectability of these objects at the new High Altitude Water Cherenkov (HAWC) 
gamma ray detector. HAWC will view the North Bubble for ~ 2 — 3 hours a day, and will map 
its spectrum at 0.1-100 TeV. For the hard primary proton spectrum required to explain five events 
at IceCube, a high significance detection at HAWC will be achieved in less than 30 days. The 
combination of results at HAWC and IceCube will substantiate the hadronic model, or constrain its 
spectral parameters. 
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Introduction. — Studying the sky at Very High Ener¬ 
gies (VHE) of ^ 100 GeV and beyond is at the frontier of 
astrophysics today. Until very recently. Gamma Ray As¬ 
tronomy has been the only avenue to probe the VHE sky, 
and has revealed the existence of powerful natural par¬ 
ticle accelerators, whose physics is still largely unknown 
[1]. It is expected that collisions of accelerated hadrons 
with the ambient medium are at least partly responsi¬ 
ble for VHE gamma ray emission from these sources. 
The same processes also produce neutrinos, and there¬ 
fore a VHE neutrino counterpart of the gamma ray flux 
is expected. In 2013, the IceCube Neutrino Observatory 
reported the detection of astrophysical neutrinos in the 
O.I-l PeV energy range mu, which opened a new av¬ 
enue of multi-messenger exploration of the VHE sky. The 
complementarity between gamma rays and neutrinos is 
strong, mainly because they both probe the sites of the 
cosmic-ray accelerators in the universe. 

Since the neutrino data are still limited by low statis¬ 
tics, a multi-messenger study of a single, specific source 
is possible only for the closest and most powerful emit¬ 
ters in the sky, like the Galactic Center and the Fermi 
Bubbles (FB). The FB have recently emerged as an ideal 
candidate for multi-messenger astronomy. These large 
globular-shaped Galactic structures, extending up to ~ 9 
kpc symmetrically out of the Galactic plane (see Fig. [^, 
were discovered in the gamma-ray data of the Fermi 
Large Area Telescope (LAT), which has since mapped 
them in details at energies ^ 0.1-400 GeV [6]. Comple¬ 
mentary observations in radio [7] and X-ray [5] confirm 
multi-wavelength emission from the FB. 

Both hadronic [91II0] and leptonic [3 HD mechanisms 
have been proposed for the origin of the FB gamma rays. 
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For hadronic models, it has been pointed out that a neu¬ 
trino counterpart should exist mm. and might be re¬ 
sponsible for a fraction of IceCube astrophysical neutri¬ 
nos [13]. In particular up to five of the IceCube astro- 
physical neutrinos of energy ^ 100 TeV-1 PeV, that are 
spatially strongly correlated with the bubble geometry 
(see Fig. |32|), are consistent with the hadronic flux 
model that reproduces the Fermi-LAT gamma-ray data 
M- Therefore, the FB could be the first object to have 
been seen in both gamma rays and neutrinos. The two 
observations, however, probe different parts of the FB 
spectrum, and therefore their connection is only indirect. 

An interesting development at the front of multi¬ 
messenger studies of the FB is expected very soon, with 
the advent of the High Altitude Water Cherenkov detec¬ 
tor (HAWC), which has started operations in Mexico this 
year (TS]. With a large field of view (~ 27r sr), excellent 
positional resolution (~ 0.1° at > 5 TeV), and strong 
sensitivity in the ^ 0.1-100 TeV range, HAWC will ob¬ 
serve the FB in great detail in a window of energy that 
partly overlaps with the regimes already probed by the 
Fermi-LAT and IceCube. It will bridge the energy gap 
between them, thus leading to a more complete mapping 
of the FB spectrum. The interplay of the Fermi-LAT, 
HAWC and IceCube data has the potential to disfavor 
or to fully establish the hadronic origin of the FB, and, 
in the latter event, constrain the spectrum of the pri¬ 
mary proton/ion flux. Such interdisciplinary study will 
be a good test bed for the development of new analysis 
techniques that could then be applied to other, less lu¬ 
minous, VHE neutrino and gamma ray sources m when 
neutrino telescopes reach the phase of high statistics data 
taking. 

The interplay between gamma ray and neutrino obser¬ 
vations of the FB is the focus of this paper. Here we 
study HAWC’s ability to detect VHE gamma rays from 
the FB in the hadronic model, and present the first dis¬ 
cussion of the potential of joint analyses of HAWC and 
IceCube data. 
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FIG. 1: Map of the Fermi Bubbles in equatorial coordinates 
(red solid contours). The blue dots show the IceCube events 
that are spatially strongly correlated with the FB with their 
sequence numbers |4] and their positional errors (blue ellipses, 
see text). Also shown are the contours of the HAWC’s field 
of view (0.6 < cos 6 < 1), for which the effective area is 
published [13, at several times of the day (magenta dashed 
contours). The shaded curves represent the inner Galaxy re¬ 
gion (corresponding to Galactic coordinates —80° < I < 80°, 
—8° < b < 8°) and the low intermediate latitude region 
(0° < I < 360°, 10° < |6| < 20°), where Fermi-L AT has 
measured diffuse Galactic gamma-ray emission m- 


Observing the FB with HAWC and IceCube. — The 
FB are extended sources in the sky subtending a total 
solid angle flpB — 0.808 sr j^. Depending on its field of 
view, location and time of the day, a detector on Earth 
will be able to observe only a fraction, /o, of this solid 
angle. Due to its privileged location near the South Pole, 
IceCube has /n — 1 [T^], meaning that it is equally sen¬ 
sitive to the FB at all times of the day. This is true for 
the High-Energy Starting Events (HESE) for which the 
vN interaction vertex is within the detector volume |3]. 

Instead, for HAWC /o is time-dependent. HAWC is 
located in the Northern hemisphere, at longitude 97.3° 
W and latitude 19.0° N. Its field of view has a half¬ 
opening angle S ~ 60°, meaning that it will be able to 
observe objects at angular distance 6 < S from its zenith 
(cos0 > cosi5 = 0.6). Fig. shows the FB in equato¬ 
rial coordinates in comparison with the time-dependent 
region observed by HAWC. We see that the north bub¬ 
ble is observable entirely at least for 2-3 hours in a day, 
while instead only a small portion of the south bubble is 
accessible to HAWC. 

For the purpose of obtaining event rates at HAWC, 
we have calculated the daily-averaged fraction of the 
solid angle (53, (/a(^1,^2)), that is subtended by the 
FB and falls in angular bins defined by the values 
cos 9 — 0.6,0.7, 0.8,0.9,1. These bins correspond to the 
intervals for which the HAWC effective are is reported 
m- For each bin, in order of increasing cos 6, we find 
{h{9i,92)) = 4.5 X 10-2, 3.5 x IO-2, 4.1 x IO-2, 1.0 x 
10-2. On average, only a few percent of flpB falls within 
a certain zenith bin. This plays the role of an efficiency 


factor in the calculation of the event rates. 

We model the expected neutrino and gamma-ray fluxes 
in a simple hadronic model, with primary proton spec¬ 
trum of the form dNp/dE oc exp{—E/Eo) and pp 
interactions with dilute gas in the FB, as described in 
details in ref. [12]. We focus on relatively hard spectra, 
with k ~ 2.2 — 2.3 and Eq ~ 3 — 30 PeV, that are re¬ 
quired by the interpretation of the five neutrino data in 
as due to the FB |14] . In Fig. the gamma ray 
and neutrino fluxes are shown for k = 2.25, Eq = 30 
PeV. Using a test, we checked that these spectra are 
a good fit of the Fermi-LAT observation of both bub¬ 
bles. In particular, ^/dof « 13/40, including penalty 
for ^2 -virhen the fit violates the upper limits (last four en¬ 
ergy bins of the Fermi-LAT data), which are treated as 
half-Gaussian. Of course, the Fermi-LAT data alone are 
also compatible with a lower-energy cutoff in the proton 
spectrum, and slightly favor Fq ^ 0(10) TeV [6]. These 
values would require a different explanation, other than 
the FB, for the neutrino data. 
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FIG. 2: The spectrum of the gamma-ray flux from the FB 
(crosses), as measured by Fermi-LAT [6]. The shaded areas 
indicate the systematic errors. The solid curve is the corre¬ 
sponding prediction of a hadronic model (parameters in the 
inset). The counterpart neutrino flux is shown too (dashed), 
and compared with the flux required by the IceGube data (the 
three dots at high energies, with error bars), in the assump¬ 
tion that FB are the sources of the five neutrino events that 
are strongly correlated with them spatially M- The neu¬ 
trino errors assume Poisson statistics, following Pearson’s 
intervals approach [21123]. The atmospheric neutrino back¬ 
ground [24], averaged over 25° — 95° zenith angle, is shown 
for comparison (steep solid line). 

Results: VHE gamma-ray and neutrino event rates. — 
As is typical of high energy cosmic-ray detectors, the 
response of HAWC and IceCube to a given flux $ (of 
gamma rays and neutrinos, respectively) can be param¬ 
eterized by an effective area, A{E,9), which depends on 
the particle arrival direction and energy. In the assump¬ 
tion that the FB are uniformly bright in gamma rays 
(and hence in neutrinos in hadronic models) [5l[6l[T2], 4) 
is zenith-independent. The number of events with pri¬ 
mary particles in a certain zenith bin, [9i , @2], and for an 





















3 


exposure time T, is given by m 

pOO 

dn / dE^E)A{E,e) 

<e<02 

POO 

^ T{fn{0i,92))nFB dE^E){AiE))e , (1) 

JEth 

where the zenith-averaged effective area (for each bin) is 
used as an approximation. Here S(<) indicates that the 
integral in the solid angle is done over the region of the 
bubbles for which the condition on the zenith angle is sat¬ 
isfied. An expression similar to eq. 0 also applies to the 
calculation of background rates in the two experiments, 
which is required to evaluate statistical significance of a 
signal from the FB. 

For neutrino event rate calculation in IceCube, we 
adopt the same method as in [14], using the averaged 
effective area as in ref. |5] for each neutrino flavor. The 
main background here is due to the flux of atmospheric 
neutrinos [23], shown in Fig.[^ The FB and atmospheric 
neutrino flux models well-reproduce the 5 cascade events 
that are strongly-correlated with the FB spatially, over a 
1000 day IceCube live time (see Fig. [^. 
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FIG. 3: Solid histogram - number of gamma ray events (in log¬ 
arithmic scale) from the FB at HAWC, in bins of the gamma 
ray energy, for 1 day exposure (left vertical axis). The corre¬ 
sponding cosmic ray background is also shown. Dashed his¬ 
togram - number of neutrino events (in linear scale) at Ice- 
Cube, in bins of neutrino energy. Here the exposure is 10^ 
days (right vertical axis), and the background is due to atmo¬ 
spheric neutrinos. The primary flux parameters are the same 
as in Fig. [^ 


For the VHE gamma-ray event rate at HAWC, both 
for signal and background, we use the effective area in 
ref. [T7] . The absorption of gamma rays due to electron- 
positron pair production with photons from starlight or 
cosmic microwave background is negligible in the HAWC 
energy range. The dominant background in HAWC is 
due to cosmic rays, mostly protons and helium nuclei, 
whose fluxes are ^4—5 orders of magnitude larger than 
the flux of the FB gamma rays. We use the hadron rejec¬ 
tion efficiency for HAWC, which is 5 x 10“^ at energies 
above 10 TeV, from ref. [2S]. Another background is 


diffuse gamma-ray emissions, which have been well mea¬ 
sured by Fermi-hAT |18| . To estimate the gamma-ray 
background in different sky regions, we adopt the spec¬ 
tra for the ^3^4^20'^ 150^'5 model (extrapolated to the 
energies of interest here) in ref. [TS] for the inner Galaxy 
region (see also ref. m for gamma-ray background in 
this region), high and low intermediate latitude region 
as shown in Fig. [^ The cosmic electron background is 
negligible. 

Figure]^ shows the signal and background daily event 
rates in HAWC from the FB region. With several sig¬ 
nal events per day, HAWC will rapidly accumulate a 
high statistics data sample. Although the signal from 
the FB is ^ 2 — 3 orders of magnitude smaller than 
the background, the high statistics will allow to estab¬ 
lish a high significance, given by the number of Gaus¬ 
sian standard deviations with respect to background as 
^ Sf/{Si -b Bi). Here Si and 

Bi are FB signal and background (total of cosmic-ray 
and diffuse gamma-ray fluxes) event rates, respectively, 
in each energy bin i. We find that, for the total of all en¬ 
ergy and zenith bins, ct > 3 (cr > 5) already after 10 days 
(35 days) of running time. After a year of operation, a 
significance of at least 5tT will be reached in each zenith 
bin separately. 

Discussion: HAWC-IceCube complementarity. — We 
have found that HAWG has an excellent potential to ob¬ 
serve the FB with high significance within a relatively 
short time scale. Depending on the parameters, this sig¬ 
nal could be consistent with the hypothesis that the bub¬ 
bles might be the (hadronic) source of the IceGube neu¬ 
trino events that spatially correlate with them. Let us 
outline below what can be learned from HAWG and Ice¬ 
Gube, in combination and individually, on the FB. 

Figure [^ elaborates on the multi-messenger connec¬ 
tion between the IceGube and HAWG, showing the re¬ 
gions of the parameter space (fc = 2.15 — 2.30 and 
Eq = 10"^ ® — 10^’® GeV) that correspond to a given sig¬ 
nal at the two detectors (significance for the HAWC and 
number of events for IceCube). For each pair {Eq, k), 
the flux normalization is determined by the low-energy 
Fermi-LAT data for the k and Eq values, similar to the 
fit shown in Fig. [^ It can be seen that, while HAWC 
can probe the entire space with high statistical signifi¬ 
cance within a year or so of operation, IceCube is in¬ 
sensitive to a neutrino flux with a spectral cutoff in the 
primary proton energy Eq ^ 1 PeV, as expected, con¬ 
sidering the higher threshold of IceCube. Therefore, the 
different combinations of possible outcomes (detection or 
exclusion) at the two observatories will be informative of 
the spectral parameters. 

IceCube has already observed 5 neutrino events that 
are correlated in position with the FB (see Fig. With 
a few more years of operation, this hint could become a 
statistically significant observation of the FB. This ob¬ 
servation would confirm the hadronic hypothesis for the 
bubbles, and indicate a primary proton population with 
a relatively hard spectrum, k ~ 2.15 — 2.25, and high 
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FIG. 4: Solid lines: contours in the space of the flux parame¬ 
ters, k and Eg, corresponding to more than 3a significance at 
HAWC (for two different running times, see legend). Dashed, 
shaded regions: isocontours of numbers of events at IceCube, 
for 1000 days running time. These numbers of events have 
been calculated at discrete points on a grid with step 0.5 
(0.05) on the horizontal (vertical) axis. 


energy cutoff, Eq ~ 3 — 30 PeV. It may provide the 
best information on the highest energy tail of the FB 
spectrum, E ^ 0.1 PeV, which is beyond the sensitiv¬ 
ity of Fermi-LAT and HAWC. Furthermore, IceCube 
will observe both bubbles, therefore testing the degree 
of symmetry between them, in a way that complements 
the lower energy information from Fermi-LAT. A neg¬ 
ative result at IceCube will be compatible with either 
a hadronic scenario with sub-PeV cutoff, or a primary 
leptonic origin of the FB. 

If IceCube observes the FB, and therefore the primary 
spectrum extends above PeV, a high statistics observa¬ 
tion in gamma rays is expected at HAWC within about 
a year of operation. This telescope is mainly sensitive to 
the North Bubble and will produce a detailed intensity 
map of it, at energies higher than those probed so far by 
the Fermi-LAT. This map will test a number of features 
seen by Fermi-LAT, like the uniform projected intensity 
of the bubbles [3E]. Due to its lower energy threshold, 
HAWC could see a signal even for a negative result at 
IceCube. In particular, thanks to its good energy resolu¬ 
tion, HAWC could reveal the presence of a spectral cutoff 
at the ^ 0.1 — 100 TeV energy window. 

Outlook. — Looking at the longer term future, a larger 
overlap between the energy windows of the two experi¬ 
ments would be desirable, especially to fully exclude the 
hadronic model in the absence of a signal at IceCube, 


accompanied by a clear detection at HAWC in the same 
energy range. This larger overlap could be achieved by 
lowering the energy threshold of the IceCube analysis. 
It has to be considered that atmospheric neutrino back¬ 
grounds at IceCube become quickly overwhelming when 
lowering the threshold, therefore the potential of such 
exercise is unclear at this time. 

Besides a direct comparison of the IceCube and HAWC 
data, the indirect influence of one experiment over the 
other will be important as well. Indeed, each data set 
will lead to improved theoretical models of the FB, which 
in turn will result in more precise predictions for the 
next generation of experiments. In particular, the ratio 
of gamma ray and neutrino fluxes may be constrained 
and compared with predictions, to help discriminate be¬ 
tween different versions of the hadronic model and possi¬ 
bly test neutrino and/or gamma ray propagation effects 
over galactic scales of distance. Similar multi-messenger 
studies for a point source or diffuse emission at the Galac¬ 
tic Center region can be done as well [26] . 

Eventually, the interdisciplinary research on the FB 
with neutrinos and gamma rays will include more par¬ 
ticipants. One of them will be KM3NeT, the multi-km^ 
neutrino detector currently planned in the Mediterranean 
sea It’s current predecessor ANTARES has already 
searched for a signal from the FB direction, resulting in a 
flux upper limit [28] • Relative to IceCube, KMSNeT will 
have a larger effective area (~ 6 km^ instrumented vol¬ 
ume) and may detect the FB in about one year of its oper¬ 
ation [27]. Finally, IceCube itself might be able to probe 
the highest energy end (> 1 PeV) of the FB gamma- 
ray spectrum using the IceTop surface detector array 
P5] . Prospects for PeV gamma-ray detection have also 
been considered for the upcoming ground-based detec¬ 
tors LHAASO (Large High Altitude Air Shower Observa¬ 
tory) [Sni and HiSCORE (Hundred*i Square-km Cosmic 
ORigin Explorer) pi] . 

We conclude that, by extending the gamma ray map¬ 
ping of the FB to 100 TeV, HAWC will be an important 
complement to neutrino searches of these objects. With 
IceCube and other neutrino observatories, it will further 
advance the new field of multi-messenger astronomy at 
very high energy. 
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